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ABSTRACT 

The role of symmetry of the inelastic light scattering amplitude, the superconducting energy gap, and the 
underlying Fermi surface manifold on the Raman spectra of unconventional superconductors is discussed in detail. 
Particular emphasis is placed on both single and bi-layer superconductors. It is found that the B\ g channel may 
be the most sensitive to doping due to the role of the Van Hove singularity. Lastly the effect of both disorder 
and spin fluctuations is considered. The theory imposes strong constraints on both the magnitude and symmetry 
of the energy gap for the bi-layer cuprates, indicating that a nearly identical energy gap of d x 2_ y 2 symmetry 
provides a best fit to the data. 
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1 Introduction 



In the recent few years Raman scattering has been proved to be a very powerful tool to study the anisotropy of 
the superconducting energy gap in high T c superconductors. The rich information provided from experiments lies 
in the polarization dependence of the obtained spectra for small light energy transfers. The recent experiments 
carried out in different high T c materials near optimal doping provide, almost identical spectra showing very 
characteristic differences between the various polarization orientations BBEI This polarization dependence, which 
is minimal in conventional superconductors, can be understood in unconventional superconductors primarily in 
terms of symmetry. 

The light scattering amplitude can be given to lowest order in the vector potential as a sum of the three terms 
which correspond to direct photon scattering via electron-hole pair production and via photon absorption and 
emission through an intermediate state. In the limit of zero momentum photon transfer this vertex can be written 
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where m is the electron mass, v stands for the intermediate state | i/k) of the electron excited out of the conduction 
band | k) with dispersion et, (^i,s an d e I,s are the incident and scattered photon energies and polarizations, 
respectively. 

According to Abrikosov and Genkin if the energy of the incident and scattered frequencies are negligible 
compared to the relevant electronic energy scale, 7 can be expressed in terms of the curvature of the conduction 
band and the incident and scattered photon polarization vectors e I,s asB 

lim n 7(k;w^ s ) = ^e;^| e s (2) 

where terms of the order of 1 — lus/luj are dropped. This effective mass approximation has been extensively used 
to calculate the Raman spectra in unconventional superconductors. 

This approximation holds only for small energy transfers and provided that the bands are not nearly degen- 
erate. Thus one might question the appropriateness of this approximation for the cuprates given that typical 
incoming laser frequencies are on the order of 2 eV - certainly on the order of the relevant electronic energy scale 
for a single band. For the case of two bands near or crossing the Fermi level, the effective mass approximation 
is even more questionable. Since in the bi-layer superconductors Bi 2:2:1:2 and Y 1:2:3 the band splitting of the 
even and odd bands is on the order of tens of meVs, the approximation misses large terms corresponding to in- 
terband transitions at relatively small energies. Detailed knowledge of the magnitude of the scattering amplitude 
in this case requires knowledge of the wave functions of the two states, which has currently not been investi- 
gated. Thereby, use of the effective mass approximation is uncontrolled in many band systems and comparisons 
to experiment intensities must be viewed with caution. 

An alternative approach is based on the experimental observation that the spectra near optimal doping for 
a wide range of cuprate materials depends only mildly on the incoming laser frequency. Since the polarization 
orientations transform as various elements of the point group of the crystal, one can use symmetry to classify the 
scattering amplitude, viz., 

7(k;wj, uj s ) = 22jl(vi, ws)$r,(k), (3) 
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where ^^(k) are either Brillouin zone (B.Z.H., orthogonal over the entire Brillouin zone) or Fermi surface (F.S.H., 
orthogonal fm the Fermi surface only) harmonics which transform according to point group transformations of 
the crystalE Representing the magnitude but not the k-dependence of both intra- and interband scattering, 
the prefactors can be approximated to be frequency independent and taken as model constants to fit absolute 
intensities. Thus we have simplified the many-band problem in terms of symmetry components which can be 
related to charge degrees of freedom on portions of the Fermi sheets. While sacrificing information pertaining to 
overall intensities, we have gained the ability to probe and compare excitations on different regions of the Fermi 
surface based solely on symmetry classifications. This can be illustrated by considering the various experimentally 
accessible polarization orientations. 

In the following we confine ourselves to tetragonal systems with in-plane lattice constant a, and for the moment 
consider uncoupled planes which can be modelled by the following band structure: 

e(k) = — 2i[cos(fc x a) + cos(fcj,a)] + At' cos(k x a) cos(k y a) — 2t"[cos(2k x a) + cos(2k y a)] — //. (4) 

We will consider coupling between the planes in the second part of the following section. Using an x, y coordinate 
system locked to the Cu02 planes, incident and scattered light polarizations aligned along x + y, x — y for example 
transform according to B\ g symmetry, and thus 

$ Blg (k) = cos^a) - cos(kya) + (5) 



where . . . are higher order B.Z.H. Likewise, e 7 ' S aligned along x, y transforms as £?2 ff : 

$B 2g (k) = sin(k x a) sm(k y a) + (6) 

The Ai g basis function is 

<&A lg (k) = ao + a2[cos(fc x a) + cos(fcj,a)] + a4COs(fe x a) cos(fc a a) + ae[cos(2k x a) + cos(2fc a a)] + . . . , (7) 

where the expansion parameters dj determined via a fitting procedure with experiment .0 The A\ g response is not 
directly accesible from experiments and must be obtained by subtracting several combinations of the response for 
various polarization orientations. 

By considering the k— dependence of the basis functions, it is clear that the B\ g part of the spectra essentially 
probes light scattering events along the k x or k y axes, Bi g probes the diagonals, and A\ g is more of an average 
over the entire Brillouin zone. It is in this manner that information^bout the momentum dependence of the 
superconducting energy gap can be usefully extracted from the data.El That information lies in two important 
aspects: the low frequency power- law behavior of the spectra and the positions of the low energy peaks. 



2 Raman Response for Unconventional Superconductors 

The spectra can be obtained via the two-particle Raman correlation function and is derivable either via Green's 
function or kinetic equation approach. In the following sections, the response is calculated for the case of a single 
band model and a superconducting bi-layer. 



2.1 Single Band 

For a single band crossing the .Fermi surface, the response can be written in the "Pair Approximation"! 
terms of the Tsuneto function A asa 
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Here E(k) 2 = e(k) 2 + A(k) 2 and we have set ks = h = 1. The light scattering cross section is obtained by taking 
the imaginary part via analytic continuation. We have neglected vertex corrections of the pairing interactions 
which while responsible for maintaining gauge invariance and producing collective modes, can be shown to have 
only a limited effect on the spectra at low energies for d— wave superconductors.^ The long-range part of Coulomb 
interaction is not included at this stage, but will be shown to be very important for the A± g response. 

Eq. (8) and approximations to it have been evaluated previously for various forms of the energy gap A(k). 
The response has been evaluated using a simpler tight binding parametrizatiap. of the energy band dispersion by 
(i) restricting the sum to the Fermi surface using a weak coupling energy gapjj rLand (ii) performing the k— sums 
directly while using a phenomenological spin fluctuation mediated energy gap.E3 These are improvements on the 
original consideration which approximated the Fermi surface to be nearly cylindrical .EJ 

In this section the T « T c Raman response for a single CuO layer is obtained by evaluating the k-sums 
directly in order to capture the charge dynamics over the entire Brillouin zoue, and in particular the role of 
the Van Hove singularities can be examined. We focus attention on La 2:l:4lLil and Y 1:2:313 as representive 
calculations. 

The B\g and Bi g responses are plotted in Fig. 1 using an energy gap of d x 2_ y 2 symmetry. In both cases, 
the Bi g response for low frequencies varies as u 3 while the B^ g varies linearly with to. This can be shown to be 



generally true for these channels for any number of higher harmonics used, and results from a consideration of 
the density of states (DOS), which varies linearly with energy due to the nodes, and the behavior of the Raman 
vertex 7 near the gap nodes. Since the B 2g vertex is finite near the nodes, the DOS determines the low frequency 
behavior. However, the gap and the vertex vanish at the same place for the B\ g channel, which in turn leads to 
the cubic behavior. This delicate interplay of vertex and energy gap can thus provide a unique determination of 
the nodal gap behavior. 

Moreover, additional information lies in the peaks of the spectra. A smooth peak is seen for B 2g and a 
double-peak is obtained for B\ g . The lower peak is due to the energy gap while the other comes from the Van 
Hove singularity. The peak of the spectra due to superconductivity in both cases is higher for the B\ g channel 
(wpeafc = 2A max ) than for B 2g , {u pea k ~ l-7A max ) since the gap maxima are located at the same part of the 
Fermi surface where the B\ g basis function is largest, and where the B 2g basis function vanishes. This is generally 
true regardless of Fermi surface shape, although the relative naaitions of the peaks can be affected slightly by 
changing the underlying manifold (see Fig 2a and 2b and RcfsJra'E]) . The shape of the spectra can also by mildly 
modified by considering higher harmonics of the energy gapl3 and/or by including final-state interactional or 
impurity scattering.E3 Apart from the last case, the power-law behaviors at low frequencies are unaffected and 
therefore their observation is a robust check on the symmetry of the energy gap. 

The Van Hove peak only shows up in the Bi g channel since the B 2g channel assigns no weight to the location 
of the Van Hove points. The Van Hove peak lies at a higher energy determined by its location off the Fermi 
surface and the value of the energy gan there. There are indications that a double-peak feature has been seen in 
underdoped Y 1:2:3 with a T c = 60KM As the Van Hove moves further away from the Fermi surface for example 
with doping, the peak shifts to higher frequencies and has a smaller residue. Inelastic scattering at higher energies 
will also act to smear out this feature and thus a single smeared and perhaps asymmetric peak would be seen if 
the Van Hove lies quite near the Fermi surface, and could in principle move the B\ g peak slightly above 2A max . 
This may be reflected in the experiments of Refill which saw a sensitivity of the B\ g peak position to doping. 

We now turn to the A\ g channel. Here the long-range Coulomb screening plays a very important role. Namely, 
those parts of the spectra which are coupled to the Coulomb fields are screened and do not appear in the limit 
of small momemtum transfer q — ► 0, as is always the case. The long-range Coulomb forces play an important 
role if there is charge transfer between unit cells at large distances. Thus only the charges produced by charge 
transfers inside the unit cell determine the measurable spectra. This transfer is between the different atoms in 
the cell or the redistribution of the electrons between different parts of the Fermi surface or between different 
Fermi surfaces. In this way, only the A\ g part of the spectra can be coupled to the Coulomb forces since for the 
case of Bi g and B 2g channels, the Coulomb corrections vanish due to symmetry. 

Including the long-range Coulomb interaction, the Ai g response can be written asBi 
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x^M = Xx7 M - , (9) 

where Xa.b is given by Eq. (8) with vertices a, b replacing 7. Writing the Ai g vertex as 7(k) = (7(k)) + 5j(k), 
where (7(k)) is the average vertex over the Brillouin zone, and inserting into Eq. (9) one finds that indeed the 
long-range Coulomb interaction completely screens the constant contribution to the vertex, or the intercell charge 
fluctuation contribution, and thus only the intracell fluctuations can scatter light in the limit of q = 0. 

In Fig. 2 we show a calculation of the unscreened and screened A\ g response for a d x 2_ y i paired superconductor 
for the Y 1:2:3 parameters. We have chosen the coefficients a, in Eq. (7) using the effective mass approximation^ 
but we note that the results do not depend dramatically on the choice of coefficient sets. By again considering 
the nature of the nodal structure of the energy gap, since the A\ g response measures a weighted average of charge 
fluctuations around the Fermi surface, it can be shown that the low energy part of the A\ g response must vary 
with energy in the same way as the DOS. 



The screened and unscreened responses at higher energies are very different. For the unscreened response, two 



peaks arise corresponding to the singular contribution of pair breaking along the axes and the Van Hove contri- 
bution as in the B\ g channel. Screening completely reorganizes the spectra to remove both peak contributions, 
replacing them with a much smoother peak at lower energies. This is due to the screening function (2nd part of 
Eq. (9)) which also contains peaks at 2A maa; due to the weighting along the k x and k y axes, and the Van Hove. 
These peaks cancel the unscreened peaks and creates a near peak lying at slightly greater than A max ^ as seen in 
Fig. 2. 

In this manner, similar symmetry considerations can be made for various types of energy gaps. The Raman 
spectra calculated for various energy gaps is summarized in Table 1 and compared to the low energy DOS N(u>). 
The dominant contribution to the Raman lineshape is due to the location and behavior of the energy gap near 
both the nodes and the maxima. Subsequently, different energy gaps produce different line shapes as summarized 
in Table 1. Here it is important to note that the cubic rise of the spectra in any channel requires that both 
the vertex and the energy gap have the same behavior near the nodes. In tetragonal systems, this requires the 
presence of a d x i_ y i energy gap. By considering small orthorhombic distortions, the A\ g and B\ g channels become 
mixed and therefore a linear rise with frequency could be obtained at low energies in a region determined by the 
amount of symmetry breaking. Since experimentally the low frequency part of the B\ g spectra rises cubically 
in the most tetragonal systems, while in Y 1:2:3 the spectra has a small linear part, this strongly suggests the 
likelihood of a d x 2_ y 2 pair state for these systems. Exploration of the low frequency part of the spectra could put 
stringent constraints on other pair-state candidates. 



2.2 Bi-Layer 

The above consideration must be reformulated for the case of two or more bands crossing the Fermi level.lll In 
particular, the role of screening is more subtle in this case. The above arguments on the role of Coulomb screening 
is based on the assumption that only one Fermi surface is relevant, jahich is very reasonable for materials with 
a single Cu02 plane in the unit cell. Recently Krantz and CardonatZI have raised the relevant question how a 
double-sheeted Fermi surface, as occurs in materials with more than one Cu02 planes in the unit cell, changes 
the above argument. In particular, the main question is whether or not in the bi-layer compound a singularity 
can occur at 2/S. max in the A% g symmetry, similar to the one which occurs in the B lg case, due to the avoidance 
of screening. 

Considering the electronic density the singular regions of the Fermi surface are near the k x and k y axes, 
respectively. The two layers double the number of these regions. In the case of B\ g symmetry the charge transfer 
changes sign going from one region to the other. In contrast, in the case of A\ g symmetry, each region has the 
same phase on a particular Fermi surface. Therefore, in the case of the single layer material, considering the 
singular part, the charge transfer is only between different cells, which is screened. In the bi-layer case, the A\ g 
symmetry allows a charge transfer between the two layers of both the same and opposite sign on each Fermi 
surface. While the addition of the two charge density distributions is still an intercell charge transfer and is thus 
screened, if the bands are non-degenerate then the difference of the two distributions can survive screening and 
could thus in principle give a singular contribution at 2A maa: . 



It is straightforward repeat the previous calculations for the two band model. The calculation proceeds by first 
considering charge density-like fluctuations on each C11O2 plane and allow for charge to be transferred between 
the planes due to an interplane coupling matrix element t± through which electrons can hop from one plane to 
the other directly or .through an intermediate state such as the chains or Y atoms. Details of the calculation 
can be found in Ref.JI^ Diagonalizing the Hamiltonian by introducing even and odd combinations of electron 
operators in the two planes, we arrive at the following compact expression for the cross section: 



„ , ( (E k ,±[7(k)±7i,2(k)]A ± (k,w) 

= g [7 (k) ± 7l , 2 (k)m ( k,.) - v __ __ _ 
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Here A± is the Tsuneto function for the bonding (anti-bonding) band, e±(k) = e(k) + ix(k), and the vertices are 
given by 

7 (k) = 71,1 (k) +72.2 (k) = i[ 7+ (k) + 7 _(k)]; 7i, 2 (k) = ±[ 7+ (k) - 7 _(k)]. (11) 

This corresponds to a diagonal contribution which allows light to be scattered by density-like fluctuations on 
either plane 1 ( 7 i,i) or plane 2 ( 7 2,2) and an off-diagonal term ( 7 i,2) which allows light scattering on both planes 
simultaneously. Since in the limit of q — > the long-range intercell fluctuations are screened, only the intracell 
fluctuations remain. Therefore only charge transfer fluctuations between the atoms in the plane and outside the 
plane (e.g., Y and the chains), and between the planes can effectively cause light scattering. Since experimentally 
the interlayer coupling is small, the resulting vertex 7 i i 2 must be smaller than the + combination, labelled as 7 . 



Rearranging terms, we can cast the result in terms of an addition of the result for each single band plus a 
mixing term: 

x» = xM + x"M + Ax», 



where 
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(12) 

(13) 
(14) 



Eqs. (12-14) shows that the total response can be considered as a sum of the single band contributions in 
addition to a mixing term which corresponds to odd combinations of fluctuations on the bands simultaneously. 
In order to discern the features of the mixing term, it is useful to write 7±(k) = 7 ^ + A 7 ±(k), where the average 
of A 7 (k) around the Fermi surface vanishes. 

Then it useful to consider the following cases: 

(1) If the layers are uncoupled, then 7 +(k) = 7 -(k) and the mixing term vanishes: Ax"(u) = 0. 

(2) If the scattering in the A\ g symmetry channel 7 is independent of k (scattering on real charge), then 
A% ~ ( 7 ° — 7 °_) 2 ~ 7 i 2- This means that the light scattering induces an interlayer charge transfer. However, it 
can be seen by the first term on the r.h.s. of Eq. (14) that in order for a divergence to occur above the peak for 
the single band result, then the energy gaps on each band must be identical. Otherwise, the divergences in the 
numerator and denominator cancel each other and a smoothly varying response is obtained. 

(3) Considering when the light induces charge transfers on different parts of the Fermi surfaces, A 7 ±(k) 
must be included and includes different symmetry channels depending on the orientation of the incident and 
scattered light polarizations. A\ g symmetry reflects the approximate tetragonal symmetry of the Cu02 planes, 
while the part of 7 showing symmetries different than A\ g (e.g., B\ gi B2 g and E g ) changes sign when symmetry 
transformations are applied. The Tsuneto function depends only on the energy gap squared, and thus all the 
sums in Eq. (14) are zero except for the A\ g component of A 7 ±(k). For tetragonal materials, this means that 
the _Bi g ,i?2g, and E g channels do not have a mixing term and therefore the total response for these channels is 
simply the additive contributions from each single band. The mixing term only contributes to the A\ g channel, 
the strength of which is roughly determined by the interlayer coupling. 

These points can be easily illustrated by considering a bi-layer superconductor with-energy gaps A±(k) = 
A±[cos(k x a) — cos(k y a)] on the two bands. For the interlayer coupling t±, we follow RefJia and define 



fj_(k) = tj_[cos(fc x a) — cos(fc y a)] 2 , 



(15) 



with tx/t = 0.2. In the absence of inter-layer coupling the chemical potential is chosen so that (n) = 0.8. Lastly, 
we consider the case where A± = 15 meV for both bands. 

The resulting spectra for the bonding and anti-bonding bands, as well as the mixing term for A\ g are shown 
in Fig. 3. Here again we have used the effective mass approximation for the expansion parameters of Eq. 7 
for illustration]!! Considering the B\ g and B^ g channels, which are simply additive, it can be inferred from 
experimental result which shows only one peak in these channels that either the energy gap is nearly identical 
both in symmetry and magnitude or alternatively that the Raman response is given predominantly by one band. 
Otherwise two distinct peaks would appear (if gaps have different symmetry and/or magnitude) and identical 
power-law behavior for each channel would be seen (if the gaps were of different symmetry but same magnitude) £3 
Since the band splitting is relatively small in Y 1:2:3 and even smaller in Bi 2:2:1:2, the second possibility in our 
opinion is unlikely. Therefore, the energy gap on each band must be similar on energy scales determined by the 
amount of inelastic scattering which would smear a double-peak feature. Again we note that the sharp nature of 
the Big peak would imply that its position and shape would be_m£>st sensitive to the location of the Van Hove 
points for each band and the maximal value of the energy gapsJijuJ 

Due to the small bi-layer splitting the mixing term for A\ g is small, as shown in Fig. (3). A very small peak at 
roughly 2A maa: (not divergent) is seen due to the odd combination of interlayer charge transfer, as suggested by 
Krantz and Cardona but is strongly surpressed. Since adding electron correlations supresses the bi-layer splitting 
even furthered the mixing term is only of minor importance and the resulting spectra can be well approximated 
as the sum of the contributions from the two individual bands. In this way, a consistent picture can emerge since 
experimentally the electronic Raman spectra do not differ substantially in the bi-layer or single-layer materials .Eraa 



3 Impurity and Spin Fluctuation Scattering 

The previous considerations were limited to the case of clean superconductors in which light scattering only 
occurs via breaking of Cooper pairs. Since the limit of q — ► was taken, as a consequence no scattering remains 
above T c due to phase space constraints. Impurity and/or spin fluctuation scattering must therefore be taken into 
account in order to attempt to model the full spectra for a wider range of temperatures and frequencies. In closing, 
this section is devoted to including such scattering processes. The following discussion is mostly qualitative in 
nature. An earlier account of this work is provided in Rcf.Jlj 

Isotropic impurity scattering can be easily incorporated using a T-matrix approach. This simplified approach 
[recent studies—have shown that impurity scattering in a strongly correlated material can have a non-trivial 
k— dependences] does not mix symmetry channels and therefore previous considerations can be applied. The 
scattering is then characterized the by the phase shift 5, c — cot (J), and by the impurity concentration n,;, 
r = Ui/nNF. Since we have seen that the shape of the Fermi surface has little effect on the low energy part of the 
cross section, we for simplicity take an infinite band and a nearly cylindrical Fermi surface [see RefB for details] . 
Here we also consider only unitary scattering c = 0, and small impuritsz_jConcentrations T/A max = 0.01. Other 
choices for the parameters are considered in a forthcoming publication.E2l While certainly a simplication of real 
systems, the salient features can be explored with some generality. 

In order to incorporate spin fluctuation scattering, we follow the work of Ref.,El who considered spin fluctuation 
scattering in a d x z_ v 2 superconductor in the context of the infrared conductivity using a 2-D Hubbard model. 
There spin fluctuation scattering was incorporated in the imaginary part of the self energy while vertex corrections 
and the real part of the self energy were neglected. We follow the same approach here and use their interpolation 
result for the frequency dependent inelastic scattering rate, which was found at low temperatures to vary as u) 
for frequencies below three times the maximum gap energy and then linearly thereafter: 

r s . / ( W ,T«T c )/3A max = 0.7( (^/ 3A ™-) 3 for a; < 3A TOOX) (lg) 

{oj/sA max otherwise 



which was obtained for T — 0.1T C , (n) = 0.85, U = 2 in a model with only nearest neighbor hopping. The 
scattering rate is approximated to be isotropic. 

The additional scattering mechanisms modify the spectra in two important regards - the low frequency behav- 
ior, caused by impurity scattering, and the high frequency behavior caused by spin fluctuation scattering. This 
is shown in Fig. (4) for the B\ g and B^ g channels. 

The impurity modifications on the density of states is reflected in the low frequency part of the spectra. Here 
the foot in the density of states causes the low frequency part of the spectra in each channel to rise linearly in 
frequency as it does in the normal state. This holds for the B\ g channel as well, where a crossover is seen at a 
frequency to* ~ y/TA max which separates linear behavior below and cubic behavior above as seen in Fig. 4. This 
is the opposite of what one would expect if the energy gap were of anisotropic— s=.wave symmetry, since there 
spectral weight would be shifted to higher frequencies due to an opening of a gapnJ'Ej The only changes at higher 
frequencies comes from a smearing of the singularities at twice the gap edge. 

Spin fluctuation scattering is relatively inoperative at low frequencies and thus does not produce any noticeable 
effects. However, the approximately linear in frequency behavior at large frequencies, as well as the value of the 
scattering rate at twice the gap edge modify the spectra at larger frequency shifts. The scattering at twice the gap 
edge effectively smears out any singularities due to the energy gap (or to the Van Hove). Further, a flat featureless 
background at larger frequencies is obtained. Similar results have been obtained by Jiang and Carbotte using a 
different model. Ell 



4 Conclusions 

We have seen that electronic Raman scattering contains a wealth of information concerning the nature of sii= 
perconductivity in unconventional superconductors which can be extracted from symmtery considerations aloneHj 
When appied to the experimental data for optimally doped cuprate materials, good agreement can be obtained 
for d x 2_ y 2 pairing in single band materials, and nearly indentical gaps of d x i_ y i symmetry for bi-layer systems. 
The analysis is based on the prominent behavior of the spectra at low frequencies and the approximate position 
of the maxima for each channel. 

However, further work is needed in order to understand the properties of the Raman spectra away from optimal 
dopingBoli3t3 It is suggested that the Bi g peak may be the most sensitive to doping due to the role of the 
Van Hove singularity. However in order to understand the effects of doping in more detail, a particular model of 
superconductivity is needed, as is a more sophisticated approach for handling the spin degrees of freedom. 
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TABLE I 



Summary of Raman response for various pair state candidates for clean, tetragonal superconductors.* 



A(k) 


N(lu -> 0) 


B\g 

x>-o) ^ 


X> - 0) 


Ai g 
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isotropic s— wave 


e(w - a) 


0(w - 2A) 2 


B(uj - 2A) 2 


0(w-2A) 2 


dx 2 — y 2 




w a 2 


to ~ 1.7 


w - 1.2 


s + id x 2_ y 2 


e(w- a.) 


9(w-2A s ) 2 


6(o;-2A s ) 2A S /A d 


6(w-2A s ) 2A S /A d 




LU 


w - 1.7 


u A 2 


w - 1.2 


(,^x 2 — y 2 ) 






w 1 /™ < 1.7 


w 1 /" 1 < 1.2 


g— wave 1 

~ cos(ft x a) cos(fcj,a) 


U) 


w 2 


w 0.6 


cj 0.6 


extended s— wave^ 
~ cos(ft x a) + cos(k y a) 


LU 


w <~ 1.7 


w 2 


w - 1.2 


experimental results 
for optimal doping 


to 


2 


w 1.5 ~ 1.8 


w 1.0 ~ 1.2 



* Here phase is undetermined and thus | d x 2 - y 2 | would yield the same as d x 2_ y 2. Moreover, orthorhombic dis- 
tortions will mix Big and A\ g channels. 

t Parameters chosen to split node at 45° into two modes at 45° ± 5° and a subsidiary maximum equal to 20 

percent of gap maxima. 

^ Parameters chosen such that (n) = 0.8. 
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Figure 1: B\ g and B2 g Raman responses for a single band plotted for parameters sets for (a) La 2:1:4 and (b) Y 
1:2:3. Here a filling (n) = 0.8 and t — 100 meV, and A (ft) = A [cos(k x a) — cos(k y a)} are used for both figures. 
The other parameters are: (a) t' /t = 0.16, t" = 0, A = 6mcV; (b) f/t = 0.2, t" /t = 0.25, A = 15 meV. 
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Figure 2: A\ g response calculated for a single band with and without Coulomb screening for the Y 1:2:3 param- 
eters. 




Figure 3: Raman response calculated for Y 1:2:3 bi-layer for various channels as indicated. The +(— ) indicates 
the bonding (anti-bonding) band, respectively, the symbol M is the mixing term [Eq. (14) which only contributes 
for the Ai g channel], and the total response is indicated by the symbol T. 




Figure 4: Raman response calculated for B\ g and B 2g channels in the absence (a) and presence (b) of impurity 
and spin fluctuation scattering. Inset shows the low frequency crossover of the B lg channel. Parameters used are 
described in text. 



